Possible existence of a band of extended states induced by inter-Landau-band mixing 
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The mixing of states with opposite chirality in quantum Hall system is shown to have delocalization 
effect. It is possible that extended states may form bands because of this mixing, as is shown through 
a numerical calculation on a two-channel network model. Based on this result, a new phase diagram 
with a narrow metallic phase separating two adjacent QH phases and/or separating a QH phase 
from the insulating phase is proposed. The data from recent non-scaling experiments is reanalyzed 
and shown that they seem to be consistent with the new phase diagram. However, due to finite-size 
effects, further study on large system size is still needed to conclude whether there are extended 
state bands in thermodynamic limit. 

PACS numbers: 73.40.Hm, 71.30.+h, 73.20.Jc 



-a 
c 

o 
o 



> 
m 

(N 

in 
o 

-I— > 



i 

C 

o 
o 



X 



I. INTRODUCTION 

According to the scaling theory of localization*, all 
electrons in a disordered two-dimensional system are lo- 
calized in the absence of a magnetic field. In the pres- 
ence of a strong magnetic field B, a series of disorder- 
broadened Landau bands (LBs) will appear, and an ex- 
tended state resides at the center of each band while 
states at other energies are localized^. The integrally 
quantized Hall plateaus (IQHP) are observed when the 
Fermi level lies in localized states, with the value of 
the Hall conductance, a xy = ne 2 /h, related to the 
number of occupied extended states(n). Many pre- 
vious ctudic j 3 ' 4 - 5 - 6 - 7 - 8 - 9 - 10 - 11 - 12 - 13 - 14 - 15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 
have been focused on so-called plateau transitions. The 
issue there is how the Hall conductance jumps from one 
quantized value to another when the Fermi level crosses 
an extended state. There are two competing proposals. 
One is the global phase diagram 5 based on the levitation 
of extended states conjectured by Khmelnitski: 2 ! and 
Laughlin* 5 ". A crucial prediction of this phase diagram 
is that an integer quantum Hall effect (IQHE) state n in 
general can only go into another IQHE states nil, and 
that a transition into an insulating state is allowed only 
from the n = 1 state. The other is so-called direct transi- 
tion phase diagram 6 in which transitions from any IQHE 
state to the insulating phase are allowed when the disor- 
der is increased at fixed B. So far, most experiments 1 *^*^ 
are consistent with the direct transition phase diagram 
although the early experiments were interpreted in terms 
of the global phase diagram. 

One important yet overlooked issue regarding IQHE is 
the nature of both plateau-plateau and plateau-insulator 
transitions. In all existing theoretical studies, these tran- 
sitions are assumed to be continuous quantum phase 



transitions. This assumption is mainly due to the early 
scaling experiments^*. The fingerprint of a continuous 
phase transition is scaling laws around the transition 
point. In the case of IQHE, it means algebraic divergence 
of longitudinal resistance slope and algebraic shrinkage of 
longitudinal resistance peak width in temperature T at 
the transition point. However, some experiments show 
that longitudinal resistance slope remains finite 2 * and 
resistance peak width remains nonzero*^ 3 - when extrap- 
olated to zero temperature. This implies a non-scaling 
behavior around a transition point, contradicting the ex- 
pectation of continuous quantum phase transitions sug- 
gested by the theories. Thus the nature of these transi- 
tions should be re-examined. 

The samples used in the non-scaling experiments 2 * are 
relatively dirty, and strong disorders should lead to a 
strong inter-Landau-band mixing. In a recent Letter—, 
we showed that the single extended state at each LB cen- 
ter broadens into a narrow band of extended states when 
the interband mixing of opposite chirality is taken into 
account. A narrow metallic phase exists between two ad- 
jacent IQHE phases and between an IQHE phase and the 
insulating phase. A plateau-plateau or plateau-insulator 
transition corresponds to two consecutive quantum phase 
transitions instead of one as suggested by existing theo- 
ries. This possibility is usually overlooked in previous 
studies where each longitudinal conductance peak is re- 
lated to a single extended state* 6 -. In this paper we shall 
present the detailed description of this study. 

The paper is organized as follows. The semiclassical 
network model for two coupled LBs is illustrated in Sec. 
II. It is shown that mixing of localized states of oppo- 
site chirality tends to delocalize a state while mixing of 
states of the same chirality does not. The similarities 
and differences between the case we consider and models 



for double layer systems and spin-degenerate systems are 
also discussed. Our approach, level-statistics technique, 
is described in Sec. III. Sec. IV has four subsections 
containing our numerical results, conclusions and discus- 
sions. Sec. IV A provides our main numerical results 
and conclusions. In Sec. IV B, possible finite size effects 
and a theoretical understanding for non-scaling behaviors 
based on standard scaling theor y 27 ; 28 are discussed. Sec. 
IV C gives further numerical results which directly show 
a narrow band of extended states in several cases, and a 
new phase diagram of quantum Hall systems is proposed. 
To further support our results, in Sec. IV D we reana- 
lyze the original data from the non-scaling experiments 21 
and show that the non-scaling behaviors in each IQHP- 
insulator transition can be attributed to two quantum 
phase transition points. Recent experiments^2i 3 A 3 i22ii& ) 
in which scaling behaviors are observed, are also dis- 
cussed and shown to be not inconsistent with early non- 
scaling experiments and our numerical results. The con- 
clusions of this paper are summarized in Sec. V. 

II. THE SEMICLASSICAL MODEL INCLUDING 
INTER-LANDAU-BAND MIXING 

According to the semiclassical theorji^i, the motion of 
an electron in a strong magnetic field and in a smooth 
random potential can be decomposed into a rapid cy- 
clotron motion and a slow drifting motion of its guiding 
center. The kinetic energy of the cyclotron motion is 
quantized by E n = (n + l/2)hio c , where uj c is the cy- 
clotron frequency and n the LB index. The trajectory 
of the drifting motion of the guiding center is thus along 
an equipotential contour of value Vq = E — E n , where 
E is the total energy of the electron. The local drifting 
velocity v(f) is determined by (in SI unit) 

v(r) = \jV{r) x B/(eB 2 ) (1) 

where S/V{r) is the local potential gradient. The equipo- 
tential contour consists of many closed loops. Neglect 
quantum tunnelling effects, each loop corresponds to tra- 
jectory of one eigenstate. The motion of electrons are 
thus confined around these loops with deviation s typi- 
cally order of the cyclotron radius l c — \Jhj{eS). 

To illustrate this semiclassical picture, let us think of 
the smooth random potential as a landscape of many 
peaks and valleys distributed randomly in the plane. 
Imagine that the landscape is filled with water up to a 
height of value Vq. The equipotential contour of value 
Vq is the boundary of land and water. According to the 
percolation theorjiS, the percolation threshold of a two- 
dimensional (2D) continuum model is p c — 1/2, where 
p c is the occupation probability of the medium (the land 
or the water). For simplicity, we suppose that the dis- 
tribution of the random potential is symmetric around 
zero. By symmetry the percolation point of both the land 
and the water is at Vb = in this case. When Vq < 0, 
the occupation probability of land is above 1/2. Thus 
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(b) 

FIG. 1: (a) Four neighboring loops in a one-band model for 
the case of Vo < 0. Dashed lines denote quantum tunnellings. 
The arrows indicate the drifting direction, (b) Loops localized 
around a valley and a peak, respectively. The arrows inside 
a loop show the directions of local potential gradient around 
the peak or valley. The arrows on a loop indicate the drifting 
direction. 



the land percolates and the water forms isolated lakes. 
These lakes are around valleys and their boundaries cor- 
respond to trajectories of localized states. In the case of 
Vb > 0, the water forms a percolating sea and the land 
becomes isolated islands around potential peaks. The 
boundary of each island is an electronic state. In short, 
semiclassical electronic states in a QH system are equipo- 
tential loops. These loops are localized around potential 
peaks for Vb > and around potential valleys for Vq < 0. 
The drifting direction of each loop is unidirectional. This 
means that they are chiral states. From Eq. ^ one can 
see that states around peaks have opposite chirality from 
states around valleys because the directions of the lo- 
cal potential gradient around a peak is opposite to that 
around a valley. If one views the plane from the direction 
opposite to the magnetic field, the drifting is clockwise 
around valleys and counter-clockwise around peaks, as 
shown in Fig^ Right at Vq = both the land and the 
water percolate, and the intersection between them is the 
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trajectory of an extended state. It means that there is 
only one extended state at Vq — for each LB. As Vo 
approaches zero from both sides, the localization length 
£ of the system tends to diverge as 

Z<x\Vo\- v (2) 

where the critical exponent v = 4/3 according to the 
classical percolation theory. Quantum effects are ignored 
in the above semiclassical argument. When two spatially 
separated loops on the same equipotential contour come 
close at saddle points of the random potential, quantum 
tunncllings should be considered. An example in the case 
of Vo < is shown in Fig^a). In the absence of inter- 
band mixing, numerical calculations have suggested that 
there is still only one extended state in each LB while the 
value of the critical exponent v is modified to be around 

7/3±a. 
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FIG. 2: (a) Two adjacent Landau bands in the case when 
the disorder broadened band width is comparable with the 
Landau gap. D(E) is the density of states. E u and E\ denote 
the centers of the two bands, (b) Schematic plot of the two 
sets of equipotential loops on the two-dimensional random 
potential for electronic states of energy E shown in (a). It 
is the projection of a three-dimensional landscape on a 2D 
vertical plane. The solid curve represents the schematic plot 
of the 2D random potential. Two dashed horizontal lines 
indicate two constant potential planes V(r) = Vo, one for the 
lower band with Vo = Vi — E — E u > and the other for the 
upper band with Vo = V u = E — E\ < 0. The ellipses denote 
the loops where the two constant potential planes intersect 
with the 2D random potential. Arrows on the loops show 
drifting directions. B is the magnetic field. 

In the case when the width of the LBs is compara- 
ble with the spacing between adjacent LBs (the Landau 
gap), and inter-Landau-band mixing should no longer be 
ignored. In order to investigate the consequences of inter- 
Landau-band mixing, we shall consider a simple system 



of two adjacent LBs. Since we are interested in interband- 
mixing of opposite chirality, we consider those states with 
energy between the lower and the upper bands whose cen- 
ters are at Ei and E u , respectively, as shown in FigEJa) . 
Thus, equipotential loops are V; = E — Ei > and 
V u = E— E u < for the lower and the upper LBs, respec- 
tively. Using the semiclassical theory described in the 
previous paragraphs, states from the upper band should 
move along equipotential loops around potential valleys 
while those from the lower band around potential peaks 
as shown in Fig^b). The loops for the upper band drift 
in clockwise direction, and those for the lower band in the 
counter-clockwise direction. These two sets of loops are 
thus spatially separated and have opposite chirality. If we 
assume that peaks and valleys of random potential form 
two coupled square lattices, the loops can be arranged 
as shown in Fig. OJa), where P and V denote peaks and 
valleys, respectively. In the absence of interband mix- 
ing, the model is reduced to two decoupled single-band 
models and all electronic states between the two LBs are 
localized. If we introduce interband mixing, the local- 
ized loops may become less localized. To see this, let 
us consider an extreme case with no tunnellings at sad- 
dle points, but with such a strong interband mixing that 
an electron will move from a loop around a valley to its 
neighboring loop around a peak and vice versa, as shown 
by B — > C in Fig[3{a). Follow the trajectory of an elec- 
tron starting at A, it will be A — ► B — > C — » D — > E ■ ■-. 
The electron is no longer confined on a closed loop, but 
delocalized! 

Before going on, we would like to make a short discus- 
sion on the relation between large disorder magnitude 
and weak magnetic field. The above semiclassical pic- 
ture for the network model is valid only when the typical 
fluctuation length of random potential, denoted by Lp, 
is much larger than the magnetic length l c — H/(eB). 
Let us take Lp as a large fixed value. Then, the effect 
of inter-band mixing is controlled mainly by two param- 
eters. One is the ratio of disorder magnitude W and 
Landau gap hu> c , which is proportional to W/B and de- 
termines whether inter-band mixing needs to be consid- 
ered. The other is the ratio of the magnetic length l c 
and the typical distance between an upper-band loop and 
its neighboring lower-band loop L p ^f- , which is propor- 
tional to W/B 3 / 2 and determines the typical strength 
of inter-band mixing in the network model. Therefore, 
the cases of strong disorder W and weak magnetic field 
B are, to some extent, equivalent in the network model 
since both of them lead to large inter-band mixing. Of 
cause, values of local potential gradients are also impor- 
tant, since they determine the detailed distribution of 
inter-band mixing in the model. 

In the one-band model, an electron can also hop from 
one loop to its neighboring loops by quantum tunnellings. 
At a first glance, this effect seems similar to that of inter- 
band mixing. However, they are fundamentally different. 
In the one-band model, electronic states for a given Vq 
have the same chirality. Thus the drifting direction of an 
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electron will be inverted when it tunnels into neighboring 
loops. This means that a strong tunnelling in a one-band 
model will induce an effective backward-scattering that 
tends to also localize electrons. We can understand this 
by considering a small part of the one-band model as 
shown in Fig^a) where all loops are moving in clock- 
wise direction. Without tunnelling, the trajectory of an 
electron starting from point A is A^B^I— > J A, 
a clockwise closed loop. With strong tunnellings, the tra- 
jectory will tend to be A^B^C~^D---^H^ A, a 
counter-clockwise closed loop. Thus, the tunnellings be- 
tween loops of the same chirality cannot delocalize states. 

Two-channel CC model can also be used to simulate 
both spin-degenerated and bilayer systems. The model 
for spin-degenerate system is the same as we use. The 
only difference is that loops for spin-up and spin-down 
states present at the same positions in real space, and 
the chiralities of loops of spin-up and spin-down states 
are the same. While the main difference between the 
model for bilayer system and ours is that the real space 
distributions of random potentials in the two layers are 
different. 

It is worthwhile to explain why we consider only those 
states between two LB centers. For states outside this 
range, both sets of loops are localized around either val- 
leys or peaks. This means that interband mixing mainly 
occurs between two loops localized around the same posi- 
tion, and this mixing will not delocalize a state. In fact, 
as explained in the previous paragraph, the mixing of 
the same chirality does not help delocalizing an electron. 
This is why we shall consider mixing between spatially 
separated states with opposite chirality. Of course, it 
does not mean that the mixing of the same chirality has 
no effect at all. As it was found in some previous works 7 , 
this kind of mixing may shift an extended state from its 
LB center. Level shifting due to mixing between states 
of the same chirality may distort the shape of the phase 
diagram, but should not alter its topology. The emer- 
gence of the bands of extended states is exclusively due 
to the mixing between states of opposite chirality. 

Now, we describe our two-channel network model in 
detail. Assume that tunnellings of two neighboring lo- 
calized states (loops) of the same band occur around 
saddle points, and interband mixing takes place only 
on the links, Fig|3fa) is topologically equivalent to the 
model shown in Fig|3Jb). FigJSfb) is the schematic illus- 
tration of our two-channel Chalker-Coddington network 
model. It is similar to the model studied in previous 
publications^^ 7 .. There are two channels on each link. 
One, denoted by a solid line, is from the lower LB around 
a potential peak. The other (dashed line) is from the up- 
per LB moving around a potential valley. The arrows in- 
dicate the drifting direction of the two sets of states. At 
each node, the tunnelling between two neighboring states 
of the same LB occurs. As shown in Fig0{a), let Z™^ 

and Z^,U be the incoming wave amplitudes of states 1 




FIG. 3: (a) Topological plot of the trajectories of the drifting 
motion of guiding centers (rhombus). The drifting motion 
around a potential peak (valley) is denoted by P (V), and 
their directions are indicated by the arrows. Dashed lines 
stand for interband mixing, and dash-dotted lines for tun- 
nelling at saddle points. The thick line (A to I) describes the 
trajectory of an electron due to a strong interband mixing, 
(b) The equivalent two-channel network model of (a). Solid 
and dashed lines on each link denote two channels from two 
LBs. Squares stand for saddle points. P, V and arrows have 
the same meaning as those in (a). 



and 2 from upper (lower) LB, respectively, and Z°^ 
and Z°^ 2 be the outgoing wave amplitudes of the two 
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states. The tunnelling is described by an S0(4) matrix 



ryOUt.l 
ryOUt,2 

\ z i / 



si 




/ 



yin.2 
ym.2 



(3) 



where the subscripts it and Z denote the upper and the 
lower bands, respectively. The elements and 
are tunnelling coefficients of an incoming wave-function 
in the upper (lower) band being scattered into outgo- 
ing channels at its left-hand and right-hand sides, re- 
spectively, s^fp, an d s^/rt are related to each other as 
„R 



i s u(i)) 2 due to the orthogonality of the ma- 



trix. Under quadratic potential barrier approximation, 
— i.e., V(x.y) = —Ux 2 + Uy 2 + V c around a saddle point, 
where U is a constant describing the strength of potential 
fluctuation and V c is the potential barrier at the point, — 
one can show that the left-hand scattering amplitude is 
given bji2& 



where e u(;) = \E + V c ■ 
electron energy, E\ = 



s u{i) = [1 + cxp(-7re M(;) )] 



-1/2 



(4) 



- (n„(Q + 1/2)E 2 ]/E 1 with E being 
f^VF^I, E 2 = ^Jk- 



1. 



and K = 



777T7TT 



+ 1. The energies of the cyclotron mo- 



tion in the two bands are (n u + 1/2)E 2 and (nj + 1/2)232, 
respectively, where U u rn are the band indices and An = 



ni 



1. The dimensionless ratio 23 r 



1 + -j^i approaches 2 from above as U or the inverse 

of uj c increases^, i.e., strong disorder regime or a weak 
magnetic field. In our calculations, we choose it to be 
2.2 since this is the regime we are interested in. For 
convenience, we choose E 2 as the energy unit and the cy- 
clotron energy of the lower band as the reference point. 
The energy regime between the two band centers is thus 

Be [0,1]. 

Inter-band mixing between two channels on a link as 
shown in FigQJb) is described by a U(2) matrix 



gout 
gout 



M 



(5) 



M 



cos v sin ( 
e 402 J \- sin 9 cos 9 J \0 e* 

' (6) 

where sin# describes the interband mixing. cf>i(i = 1 ~ 4) 
are random Aharonov-Bohm phases accumulated along 
propagation paths. In our calculations, we shall as- 
sume that they are uniformly distributed in [0, 27r)24. 
In the following discussion, a parameter J, defined as 
\J J/(l + J) = sin 9, is used to characterize the mixing 
strength. J will take the same value for all links in our 
calculations. We hope that this simplification will not 
affect the physics. 
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FIG. 4: (a) A node with four incoming and four outgoing 
channels. Z^f, is the wavefunction amplitude of the i" 1 in- 
coming wave from the upper (lower) LB. is that of out- 
going wavefunction amplitude, (b) A link with two channels. 

is the incoming (outgoing) wavefunction amplitude 
of the upper (lower) LB. 



III. THE APPLICATION OF 
LEVEL-STATISTICS TECHNIQUE ON THE 
NETWORK MODEL 

Electron localization length is often obtained from the 
transfer matrix method. For a two-dimensional system, 
however, it is well known that this quantity alone does 
not provide conclusive answers to questions related to 
the metal-insulator transition (MIT)^. On the other 
hand, level-statistics analysis has been used in study- 
ing MIT^S^i. Level-statistics analysis is based on ran- 
dom matrix theory (RMT) 4 — . The basic idea is that the 
localization property of an electronic state can be de- 
termined by the statistical distribution function P(s) of 
the spacing s of two neighboring levels. For localized 
states, the distribution is Poissonian Ppe{s) — exp{—s), 
called 'Poissonian ensemble (PE)'. In the case of extended 
states, the nearest neighbor level spacing distribution has 
the following formic 



P{s) = Ci{P)sPexp[-C 2 {P)i 



(7) 



where C\{j3) and C 2 (j3) are normalization factors deter- 
mined by J P(s)ds — 1 and J sP(s)ds = 1. The pa- 
rameter [3 is determined by the dynamical symmetry of 
the system. The case of j3 = 1 is for systems with time- 
reversal symmetry and an integer total angular momen- 
tum and is referred as 'Gaussian orthogonal ensemble'. 
Systems with time-reversal symmetry and a half-integer 
total angular momentum belong to the case of j3 = 4, 
called 'Gaussian symplectic ensemble'. For systems with- 
out time-reversal symmetry (3 = 2, and it is called 'Gaus- 
sian unitary ensemble (GUE)'. A fundamental difference 
between level statistics property of extended states and 
localized states is that P(s) at s = is zero for extended 
states and one for localized states. The physical reason 
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of this difference is the so-called 'level repulsion' of ex- 
tended states. Two extended states with the same 'bare 
energy' will overlap in real space(since they are extended 
in real space) and form two new extended states with dif- 
ferent energies. While localized states can have the same 
energy staying in different regions of real space. This 
approach is appropriate for the network model of quan- 
tum Hall systems because localized states in the model 
are loops at different regions of real space while extended 
states in the model are formed by quantum percolation 
of such localized loops. 

We shall follow the approach proposed by Klcssc and 
Metzler 43 . A quantum state of a network model can be 
expressed by a vector whose components are electronic 
wave-function amplitudes on the links. In our case, the 
vector can be written as $ = ({Z^Zf}), where Z l u and 
Z\ are the electron wave-function amplitudes of the up- 
per band (it) and the lower band (0 on the i-th link, 
respectively. As shown by Fertig 4 ^, the network model 
can be described by an evolution operator U(E), an in- 
dependent matrix determined by the scattering proper- 
ties of nodes and links in the model. (As an example, 
the evolution operator of a two-channel network of size 
L = 2 with periodic boundaries on both directions is 
constructed explicitly in the Appendix.) In general, the 
eigenvalue equation of the evolution operator is 



step can be written in the following form 



*v(t- 





Un- 



ify- 

. v 6 



• H 



(9) 



U(E)$ a (E) = e^ E ^ a (E), 



(8) 



where is the (2L 2 ) x {21?) zero matrix. Uv^h de- 
scribes how wavefunction on vertical links evolves into 
that on the horizontal links. Similarly, Uh^v describes 
that from horizontal to vertical links. For the detail 
derivation, we refer readers to the example shown in the 
Appendix. The evolution equation in two time steps is 
given as 

$ H (< + 2) = U v ^ffUff^ v $ H (t) (10) 
$v(t + 2) = U H ^vUv^H$v(t) (11) 

Therefore, the evolution matrix in two time steps is 
block-diagonal and the two blocks have essentially the 
same statistical property. We thus need only consider 
one of them. 

We study the model for 1=8, 12, 16, 20 and 24. The 
calculation procedure is as follows. Take a realization of 
the random phases, construct the evolution matrix and 
obtain the quasi-energies {tOi}. Put them in descending 
order and calculate the level spacings Si = {uji — Ui—i)/6, 
where S is the average of Sj. Repeat this procedure for 
sufficient times so that more than 5 x 10 4 level spacings 
are collected for a given E and J. The level-spacing 
distribution function P(s) is thus obtained numerically. 



where a is the eigenstate index of U. The true eigenen- 
ergies {E n } of the system are those energies at which 
uj a (E) is an integer multiple of 2n. It has been shown 
by Klesse and Metzler— that the set of quasi-energies 
{u> a (E)} corresponds to the excitation spectrum of the 
stationary state with energy E. Therefore, the statistics 
property of the set of quasi-energies {uj a (E)} at E is the 
same as that of true eigenenergies {E n } around E, and 
the localization property of an electronic state with an 
energy E can be obtained by the quasi-energies. The 
advantage of this approach is that all the quasi-energies 
can be used in the analysis so that better statistics can 
be obtained. 

Chalker and CoddingtonM showed numerically that 
an open boundary condition along one direction creates 
extended edge states along the other direction. In or- 
der to get rid of the edge states, we employ a periodic 
boundary along both directions in our calculation. For a 
two-channel network model of L x L nodes with periodic 
boundaries along both directions, there are 4L 2 compo- 
nents in <I>. U is thus a (4L 2 ) x (4L 2 ) matrix. However, 
there is a special property of the network modelfe the 
nodes scatter electrons only from vertical channels into 
horizontal channels and vice versa. If one separates $ 
into the set of wavefunction amplitudes on the horizon- 
tal links and the set of wavefunction amplitudes on 
the vertical links the evolution equation in one time 



IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. Analysis of the level-spacing spectrum 

We shall analyze the numerical results of the level- 
spacing distribution function P(s) in order to provide 
evidences for the existence of extended state bands in 
our model. Due to the chirality nature of the drifting 
motion, time-reversal symmetry is absent from our semi- 
classical network model. Then, according to the RMT— , 
P(s) should be the GUE distribution, Pgue{s) = 
32ir~ 2 s 2 exp(—4:S 2 /n), for extended states, and the PE 
distribution Ppe(s) for localized states. Since the over- 
all shape of the GUE distribution is quite different from 
that of the PE distribution, one may use P(s) to dis- 
tinguish an extended state from a localized state. Fig0 
is P G ue{s) and P(s) for (E = 0, J = 0.1) (a), (E = 
0.02, J = 0.1) (b), and (E = 0.5, J = 1.5) (c) with 
L = 8, 12, 16, 20, 24. Figgis P(s) for (E = 0.0, J = 0.7) 
(a), (E = 0.02, J = 0.7) (b) and (E = 0.5, J = 0.5) (c). 

The overall shape of these curves has some common 
features. All curves have a vanishing value when s goes to 
zero. At small s they increase with s and reach a peak at 
some intermediate s. Then they decrease monotonically 
to zero with increasing s. These features are the same 
as those for Pq UE (s)^. Thus they all look like to be 
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FIG. 5: P(s) vs. s for L = 8, 24 and Pous(s). (a) £ = and 
J = 0.1; (b) E = 0.02 and J = 0.1; (c) E = 0.5 and J = 1.5. 



Pgue(s) at first glance. This raises the question of how 
to distinguish numerically extended states from localized 
states. As a simple way, it is natural to expect that P(s) 
for an extended state approaches Pgue{s) while that for 
a localized state should deviate from Pgue{s) as L in- 
creases. Indeed, as L increases, curves in each sub-figure 
of Fig O approach Pgue(s) while those in FigHjJshow the 
opposite tendency. Thus we can use this different ten- 
dency of P(s) to distinguish extended states from local- 
ized states. We shall show quantitatively that such op- 
posite tendency for extended states and localized states 
exists also in several other characteristic quantities. 

Let us first consider a characteristic quantity Iq de- 
fined by Iq = J s 2 P(s)ds/2. It is commonly used to 
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FIG. 6: P(s) vs. s for L = 8, 24 and Pgue(s). (a) E = and 
J = 0.7; (b) E = 0.02 and J = 0.7; (c) E = 0.5 and J = 0.5. 



characterize the overall shape of P(s) and to examine 
the localization property^. It is well-known that Iq = 1 
for localized states while To < 1 for extended states^. 
Thus, the following simple criteria is employed: A state 
is localized if its Iq increases and approaches 1 as I in- 
creases. Otherwise, it is extended. Curves in Fig0 arc 
Jo vs. mixing strength J for E = (a); 0.02 (b); and 0.5 
(c) for L = 8, 12, 16, 20, 24. Figl^b) shows that the state 
of E = 0.02 is localized at zero mixing and extended at 
small J. Then it is localized again after J passes a par- 
ticular J c where Iq of different L cross. For the state of 
E = at the lower band center shown in FigEJa), it is 
extended at zero mixing. Then, it shows the same feature 
as the state of E = 0.02 at small and large J. Fig[7Jc) 
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FIG. 7: Jo vs. J for L 
E = 0.02; (c) E = 0.5. 



8,12,16,20,24, (a) E = 0; (b) 



FIG. 8: I P (s) vs. s for £ 
Pgc/e(s). (a) E = and J = 
(c) £ = 0.5 and J = 1.5. 



= 8, 12, 16, 20, 24 and that for 
0.1; (b) E = 0.02 and J = 0.1; 



shows that state of E = 0.5 is always localized at small 
J and extended only for large J(> 1) where all curves of 
different system sizes tend to merge together. 

A fundamental difference between P(s) for a localized 
state and that for an extended states is its behavior at 
small s. As s goes to zero, P{s) vanishes for extended 
states due to level-repulsion while it approaches 1 for lo- 
calized states due to level-aggregation^. Thus we need 
to consider the behavior of P(s) at small s for further 
test of the results in the last paragraph. It is convenient 
to consider a function of integrated level-spacing distri- 
bution at small s defined by Ip(s) — J* P(s')ds' . Ip(s) 
is the fraction of level-spacing smaller than s. Although 
P(s) in most cases of our numerical results is close to 



the GUE distribution, level-repulsion of extended states 
and level-aggregation of localized states should still be 
expected at small s. This leads to the following cri- 
teria: Ip(s) at small s should increase (decrease) with 
L for localized (extended) states. Thus the behavior of 
Ip(s) at small s can serve as another way of distinguish- 
ing extended states from localized ones. FiglHl shows 
Ip(s) for [E = 0,J = 0.1) (a), (E = 0.02, J = 0.1) 
(b) and (E = 0.5, J = 1.5) (c) for L = 8,12,16,20,24 
and comparison with Ip(s) of Pgue(s)- FigEj] is for 
(E = 0,J = 0.7) (a), (E = 0.02, J = 0.7) (b) and 
(E = 0.5, J = 0.5) (c). One can see clearly that states 
in FigJH] show the feature of extended states while states 
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FIG. 9: Ip(s) vs. s for L = 8,12,16,20,24 and that for FIG. 10: 7p(s = 0.5, J) vs. J for £ = 8,12,16,20,24, (a) 
Pgue(s), (a) E = and J = 0.7; (b) E = 0.02 and J = 0.7; J5 = 0; (b) E = 0.02; (c) E = 0.5. 
(c) E = 0.5 and J = 0.5. 



in Fig|5|are localized. In order to examine an electronic 
state of fixed energy in the whole range of mixing, we 
consider Ip(s = 0.5), the fraction of the level-spacings 
less than 0.5. We plot the results of Ip(s = 0.5) vs. J 
at E = 0,0.02 and 0.5 for L=8, 12, 16, 20, 24 in FigHDJ 
Similar to the criteria for Ip(s), we use the following 
ones. If Ip(0.5) of a state increases with L, the state is 
localized. Otherwise, it is extended. According to this 
criteria, curves in Fig llOl give essentially the same results 
as those obtained from Iq in Fig[7| 

Let us now turn to the region of large s. Since Pgue (s) 
decays faster than Ppe(s) at large s, the behavior of P(s) 
in this region can also be used to differentiate extended 



and localized states. In this region it is convenient to con- 
sider another quantity defined by F(s) = f P(s)ds = 
1 — Ip(s). The meaning of F(s) is the integrated fraction 
of level-spacings larger than s. Since F(s) of Pgue( s ) is 
less than that of Ppe(s) at large s, we may expect that 
F(s) at larger s decreases (increases) with L for extended 
(localized) states. Fig^Jis F(s) for Pgue(s) and that 
for {E = 0, J = 0.1) (a), (E = 0.02, J = 0.1) (b), and 
(E = 0.5, J = 1.5) (c) with L = 8, 12, 16, 20, 24. Figgis 
F(s) vs. s for (E = 0, J = 0.7) (a), (E = 0.02, J = 0.7) 
(b), and (E = 0.5, J = 0.5) (c). In view of Fig and 
Fig|51 it is clear that the results of F(s) coincide with 
those of Ip(s) concerning the localization property. We 
also calculate F(s — 2) for fixed energy states in the 
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FIG. 11: F(s) vs. s for L = 8,12,16,20,24, and that for 
Pgue(s). (a) E = and J = 0.1; (b) E = 0.02 and J = 0.1; 
(c) E = 0.5 and J = 1.5. 



FIG. 12: F(s) vs. s for L = 8, 12, 16, 20, 24 and that for 
Pgue(s). (a) E = and J = 0.7; (b) £ = 0.02 and J = 0.7; 
(c) E = 0.5 and J = 0.5. 



whole range of mixing. As shown above, the same cri- 
teria as that for Iq and Ip(s = 0.5) can be employed. 
The curves of F(s = 2) vs. J are plotted in Fig 1 131 for 
E = (a), E = 0.02 (b) and £? = 0.5 (c). One can see 
that they are consistent with the results of Jo (FigEJ and 
I P (s = 0.5) (FigJTDJl. 



B. Discussion of finite-size effect 

In this subsection, we shall consider possible influence 
of finite-size effect on our numerical results. It is known 
that localization lengths for 2D models can exceed 3 x 10 4 



lattice spacings, a thousand times larger than the max- 
imum lattice size L — 24 in our numerical calculation. 
Therefore, one should worry about finite-size effect, and 
question the unsuitability of our simple criteria for lo- 
calization property. In the following discussion, we shall 
only examine the results of the quantity Iq- Essentially 
the same discussions can be made for other quantities 
such as Ip{s) and F(s). 

Let us first consider the results for E — and E = 0.02 
in Fig[7| The two cases are quite similar. Curves for dif- 
ferent L cross at a single point J c . As L increases, Iq 
decreases and approaches the value for extended states 
when J < J Cl while it increases and approaches the value 
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FIG. 13: F(s = 2, J) vs. J for L = 8, 12, 16, 20, 24, (a) E = 0; 
(b) £ = 0.02; (c) E = 0.5. 



for localized states when J > J c . A straightforward way 
of interpreting this behavior is that the state exhibits 
a transition from an extended state in J < J c to a lo- 
calized one in J > J c . The correlation length diverges 
at critical point J — J c and decreases sharply when J is 
slightly away from J c . (In a metallic phase the correlation 
length is small while the localization length is divergent.) 
Fluctuations at length scale of order of the correlation 
length cause Iq to deviate from its thermodynamic-limit 
value, but Jo approaches its thermodynamic-limit values 
for both J < J c and J > J c as lattice size increases. This 
is a natural explanation of the results. The only finite- 
size effect in this explanation is rounding behavior in a 
region close to the critical point J = J c . This region is 



normally very narrow because of the sharp drop of the 
correlation length near J = J c . 

Another possible interpretation is to assume that 
states are always localized in the whole region except at 
J = J c . In this case, the localization length equals to the 
correlation length. Then, in order to explain the above 
behavior, we have to assume the following features of the 
localization length: (a) In the region J < J c , the ratio of 
localization length £(L) and system size L should increase 
with L for small L, leading to a false metallic behavior, 
while this ratio decreases with L for large L, a behavior 
for the localized state; (b) in the region J > J c , £(L)/L 
should always decrease with L. In principle, one cannot 
rule out this possibility without doing calculations for 
large lattice sizes. However, (a) is too strange to justify. 
Furthermore, there is no reason to believe why (a) occurs 
in region J < J c , but not in region J > J c . 

For E — 0.5 in Fig[7| localized behavior is clearly seen 
for small and intermediate J while the curves of different 
system size tend to merge at large J. Here the finite- 
size effect should be considered seriously. There are two 
ways to explain the merging behavior. One is that a line 
of critical points exists for large J where the correlation 
length is always divergent. The other is that the correla- 
tion length at the thermodynamic limit is very large yet 
finite and the merging behavior is just a finite-size effect. 
Unlike the case of E = and E = 0.02, both explana- 
tions in this case are reasonable. The only way to make 
an unambiguous conclusion is to do calculations for large 
sizes. 

However, we can propose a physical picture for the ex- 
istence of new extended states at E ~ 0.5 in the case 
of strong interband mixing. Assuming that the intra- 
band tunnelling at nodes are negligibly weak for states 
of E ~ 0.5, we saw already from Fig. Ola) that the 
maximum interband mixing (sin# = 1) delocalizes the 
state, which is localized at zero interband mixing. If one 
views p = sin 2 9 as connection probability of two neigh- 
boring loops of opposite chirality, our two-channel model 
without intra-band tunnellings at nodes is analogous to 
a bond-percolation problem. It is well-known that a per- 
colation cluster exists at p > p c = 1/2 or J > J c = 1 
for a square latticed. Therefore, an extended state is 
formed by strong mixing. One hopes that the intra-band 
tunnellings at nodes will only modify the threshold value 
of the mixing strength. If this picture is correct, the 
finite-size effect only affects our efforts to determine the 
accurate value of J c yet it does not influence the existence 
of such a critical point. 

It should be noted that all above discussions are based 
on the single parameter scaling argument. Suppose that 
the region of extended states in thermodynamic limit is 
vanishing instead of remaining finite, our above numeri- 
cal results can also be explained by introducing irrelevant 
length scales and considering their corrections to scaling, 
as has been pointed out by Pruiskenl and Hucksteir. 28 . 
Thus both our picture of finite region of extended states 
and the idea of corrections to scaling are alternative ex- 
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planations for non-scaling behaviors. Therefore, calcu- 
lations for larger sizes of system are necessary to make 
unambiguous conclusion for the width of extended state 
region in thermodynamic limit. In a recent theoretical 
study, Pruisken et. ali^£ have developed a microscopic 
theory based on nonlinear a model to explain the non- 
scaling behaviors within the assumption of a single crit- 
ical point. However, in their consideration of the case of 
long-range correlated disorder which corresponds to the 
network model, interband mixing, the physical reason for 
the existence of finite region of extended states in our pic- 
ture, is neglected. Therefore, their study does not rule 
out the possibility of finite region of extended states. 
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C. Discussion of localization property 

In the past two subsections, we analyzed the global 
shape of P(s) and its behavior at small and large s by 
considering Iq, Ip(s) and F(s), respectively, and discuss 
possible influence of finite-size effect. Analysis of all these 
quantities leads to essentially the same conclusion con- 
cerning the localization property, as follows. For zero 
interband mixing, only states at the two LB centers are 
extended. In the presence of interband mixing, new ex- 
tended states emerge. States near the LB centers, — i.e., 
E ~ 0, — are delocalized by weak interband mixing and 
localized by strong mixing, with a transition point at 
some intermediate mixing J c . For states near the region 
between two LBs, — i.e., E ~ 0.5, — they are localized at 
both weak and intermediate mixing and delocalized by 
strong mixing. 

In order to show explicitly the existence of a narrow 
band of extended states and its evolution with increasing 
mixing, curves of Iq vs. E are plotted for three values of J 
in Fig^l A band of extended states is formed around the 
LB center EE [0, 0.1] for J = 0.1. When J is increased to 
an intermediate value 0.7, this band of extended states is 
lifted up to E e [0.8, 1.6]. For strong mixing, it is further 
shifted to E e [0.4,0.5]. Thus the band of extended 
states in the lower LB tends to float up in energy while 
that of the upper tends to dive down in energy as mixing 
strength increases. The two bands should finally meet at 
the middle energy region in the case of strong mixing. 

The above results are restricted to the case of two LBs. 
However, there are infinite number of LBs in a realistic 
system. In order to conjecture the situation when infinite 
number of LBs is incorporated in our result, we take into 
account the float-up-merge picture proposed by Sheng 
et. alii£. We shall expect that a narrow extended band 
appears in each LB center for weak interband mixing. 
Increasing mixing, i.e., increasing disorders or decreasing 
magnetic field, the extended band in the lowest LB floats 
up and finally merges with that in the second lowest LB. 
Then, this extended band will further shift up and merge 
into that in the third lowest LB, and so on so forth. 

To express our numerical results in the plane of en- 
ergy and interband mixing, a topological phase diagram 
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FIG. 14: J vs. E for L = 8,12,16,20, (a) J = 0.1; (b) 
J = 0.7; (c) J = 1.5. 



shown in Fig ll5f a'l is obtained. In the absence of inter- 
band mixing, only the singular energy level at each LB 
center is extended. In the presence of interband mix- 
ing of opposite chirality, there are two regimes. At weak 
mixing, each of the extended states broadens into a nar- 
row band of extended states near the LB centers. With 
increased mixing, the extended states in the lowest LB 
shift from the LB center(see Fig ll4(l . These extended 
states will eventually merge with those from the higher 
LBs. This shifting of bands of extended states is similar 
to the shifting of single extended states at LB centers ob- 
served in previous studies 5 where emerging of extended 
bands is missing. At strong mixing, a band of extended 
states exists between neighboring LBs where all states 
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are localized without the mixing. 

Let us look at the consequences of the above results. 
For weakly disordered systems in IQHE regime, the Lan- 
dau gap is larger than the LB bandwidth. Thus there is 
no overlap between adjacent LBs. According to the semi- 
classical picture, electronic states between the two adja- 
cent LBs should be from either the upper or the lower 
bands with the same chirality in this case. It means that 
no interband mixing occurs and there is only one ex- 
tended state in each LB. This may explain why scaling 
behaviors were observed for plateau transitions in early 
experiments on clean samples. Interband mixing occurs 
when the Landau gap is less than the LB bandwidth. Sys- 
tems of relatively strong disorders in IQHE regime should 
correspond to this case. As the single extended state at 
each LB center broadens into a narrow extended band, 
a narrow metallic phase emerges between two neighbor- 
ing IQHE phases. Thus each plateau transition contains 
two consecutive quantum phase transitions for strongly 
disordered systems. The bands of extended states will 
merge together in strong mixing. This strong mixing 
regime corresponds to the case when the Landau gap is 
much smaller than the bandwidth. Since the Landau 
gap is proportional to the magnetic field, the disordered 
system should always enter the strong mixing regime be- 
fore it reaches the weak field insulating phase, regardless 
of how weak the disorders are. In terms of QH plateau 
transitions, a direct transition occurs because a narrow 
metallic phase exists between two QH phases in a weak 
field. Thus, we propose that a direct transition from an 
IQHE phase to the insulating phase at weak field is real- 
ized by passing through a metallic phase, and it should 
hold for both weak and strong disordered systems. 




the weak localization theoryi, no extended state exists 
at this point. Differing from existing theories, there ex- 
ists a narrow metallic phase between two adjacent IQHE 
phases and between an IQHE phase and an insulating 
phase. 

D. Comparison with previous studies 

In this subsection, we shall compare our results with 
previous studies. We shall show that the new phase 
diagram in Fig ll5l is consistent with the non-scaling 
experiments— where samples are relatively dirty, and 
interband mixing is strong, corresponding to a process 
along line a in Fig^Jb). The system undergoes two 
quantum phase transitions each time when it moves from 
the QH insulating phase to IQHE phase of n = 1 and 
back to the weak field insulating phase as the magnetic 
field decreases. To verify this claim, we analyzed the 
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FIG. 16: The original experimental data of 

\D.[Rxx{f,T)/R* x (0.647, 40mK)\ in Ref. 18. / is the 
filling factor of LBs and T is the temperature. 

original experimental data in Ref. 18 according to the 
assumption of two quantum phase transition points. The 
experiment data of the logarithm of the longitudinal re- 
sistance ln[R xx (f,T)] are shown in FiglTol where / is the 
filling factor of LBs and T is the temperature. According 
to the theory of continuous transitions, one should obtain 



FIG. 15: (a) Topological phase diagram of electron localiza- 
tion in E — J plane. The shadowed regime is for extended 
states (metallic phase) . (b) Topological QH phase diagram in 
W — B plane. W stands for the disorder strength, and B for 
the magnetic field. The shadowed regime is for the metallic 
phase. The area indicated by the symbol n is the n-plateau 
IQHE phase. The rest area is the insulating phase. 

Plotting above results in the plane of disorder and the 
magnetic field, we obtain a new topological QH phase di- 
agram as shown in Fig ll5f . This is similar to the em- 
pirical diagram obtained experimentally in Ref. 15. The 
origin (W = 0, B = 0) is a singular point. According to 



ln[R xx (u,T)] = F 1 (S 1 (f)/T) (12) 
with Sitf) ~ (fd - f) ZlUl for the region of / < f cl while 

ln[R xx (u,T)]=F 2 (S 2 (f)/T) (13) 

with S 2 (f) ~ (/ ~ fc2) Z2U2 for the region of / > f c2 . 
Previous theories predict one single critical point, — i.e., 
fci = fc2 and Z\V\ = z 2 v 2 . But our results suggest two 
distinct critical points. By standard scaling analysis, two 
good scaling behaviors are obtained for two close critical 
filling factors of f c \ = 0.6453 and f c2 = 0.6477 as shown 
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FIG. 17: The fitting result of two critical points at the left and 
the right side. The two straight lines show coincidence with 
the scaling law. The critical filling factors are f c i ~ 0.6453 
and f C 2 = 0.6477. The two critical exponents are equal to the 
value zv = 2.33 ± 0.01. 
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FIG. 18: The best fitting result of one single critical point at 
the left and right side. The critical filling factor is f c = 0.646. 
The two straight lines illustrate systematic deviations from 
the scaling law at regions close to the critical point. The 
average values of the two critical exponents are z\ v\ = 2.58 ± 
0.02 and z 2 v 2 = 2.60 ± 0.02, respectively. 



in Figll7l The critical exponents in both the left side 
and the right side of the transition region are equal to 
the value zv — 2.33 ± 0.01. On the other hand, the fit 
for one single critical point fails. Fig^3] shows the result 
of a single critical point at v c — 0.646. It is the best 
fitting result for a single critical point if we require that 
the two critical exponents are approximately equal and 
the scaling law is optimally obeyed. The two critical 
exponents are z\Vi = 2.58 ± 0.02 and 2:2^2 = 2.60 ± 
0.02, deviating from the theoretical results zv ~ 2.33. 
One can also see clearly systematic deviations from the 
scaling law in the region close to the critical point at both 
sides in Fig^| This implies that the transition process 
is governed by two separated critical points instead of 
one. The regime between the two critical points should 
correspond to the metallic phase. 



It is worth noting that there is another puzzle in the 
non-scaling experiment which may be solved by our two- 
critical-point picture. As an example, we consider the 
experimental data for the transition between the QH 
insulating phase and the n — 1 IQHE phase. It was 
shownSi that the logarithm of the longitudinal resistance 
ln[R xx (f, T)] can be fitted by a linear function of the 
filling factor / (see Fig ll6fl 

ln[R xx (f,T)]=ln[R xx (f c ,T)]-(f-f c )/(a + pT)] (14) 

where a and (3 are positive constants, f c is the filling fac- 
tor where curves of different temperature T cross approx- 
imately. Since a is non-zero^i, it leads to the conclusion 
that R X x(f,T) at the limit of T = remains finite for 
every f. This is puzzling because it is inconsistent with 
the theoretical requirement that R XX (T = 0) = 00 in the 
QH insulating phase, i.e., f < f c , and R XX (T = 0) = 
in the n = 1 IQHE phase, i.e., / > f c . This puzzle may 
be solved as follows. Combine the linear relationship be- 
tween ln[R xx (f, T)] and / for fixed T with our picture of 
two critical points f c ± < / C 2, we expect 

ln[R xx (f,T)]=ln[R xx (f cl) T)}-(f-f cl )/(A 1 T z ») (15) 

in the QH insulating phase, i.e., / < f c \, while 

ln[R xx (f,T))=ln[R xx (f c2 ,T)]-(f-f e2 )/(A 2 T z ») (16) 

in the n=l IQHE phase, i.e., / > f c2 , where A\ and A 2 
are positive constants, and z and v are critical exponents. 
It is clear that both R X x{f, T = 0) = 00 in / < f c i and 
Rxx(f,T = 0) = in / > f c2 are recovered. While a 
finite value of R xx (/, T = 0) in the region f c i < / < 
f c2 is consistent with our prediction of a metallic phase 
between the two critical points. 

Scaling behaviors of plateau-plateau transitions are 
observed in recent experiment922^2i2ii^2iS. At a first 
glance, it seems that these results conflict with both 
non-scaling experiments and our numerical results. How- 
ever, there are two important differences between re- 
cent scaling experiments and non-scaling experiments. 
One is that some scaling experiments22^2i2i were done 
deep inside QH plateau, or far from the plateau transi- 
tion point, while the non-scaling behavior was obtained 
by using data very close to the transition points. In 
the non-scaling experiments 2 ' 1 , it is known that data 
not too close to a transition point follow a scaling 
law. The other is that the samples used in all scaling 
experiments22i22i2ii22iS are clean with very high mobil- 
ity while non-scaling behavior was observed in relative 
dirty samples2i*22i2£. In fact, the mobility in recent scal- 
ing experiments is at least one order of magnitude larger 
than that in early non-scaling experiments. This means 
that the scaling and non-scaling experiments correspond 
to regions of vanishing and relatively strong inter-band 
mixing, respectively. Thus, there is no real conflict be- 
tween the recent scaling experiments and those early non- 
scaling experiments. Since our model is valid only when 
inter-band mixing is considerable, there is also no real 
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conflict between recent scaling experiments and our nu- 
merical results. 

Two-channel CC model has been used to simulate two 
degenerated or nearly-degenerated spin-resolved Landau 
subbands with strong interband mixing by Wang et ali. 
They found two distinct critical points which were re- 
lated to mixing- induced repulsion 7 . For the nearly- 
degenerated case, they did not consider the states be- 
tween the two LB centers. For the degenerated case, 
their study could not discern whether the states between 
the two critical points are extended or localized. Accord- 
ing to our results, a band of extended states is formed 
around the degenerated LB center for the degenerated 
case at strong mixing, while for the nearly-degenerated 
case electronic states between the two LB centers are 
delocalized by strong mixing. Thus our results are con- 
sistent with their results. 

One should also notice that two types of metallic 
phases have been studied extensively in the QH system. 
One is the composite Fermion state at the half-filling 
in the lowest Landau level (LL) and the other is the 
stripe state at the half-filled higher LLs. These states 
are formed by the Coulomb interaction effect in the high 
mobility samples. They are different from our metal- 
lic phase caused by level mixing. Although we have not 
considered the electron-electron interactions in our study, 
there is no reason why the delocalization effect of level 
mixing will be diminished by the Coulomb interaction. 
Of course, the interaction could lead to a level shift, thus 
it may modify the band width. 
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CONCLUSIONS 



In conclusion, we find by numerical calculations within 
the network model that it is possible that the single ex- 
tended state at each LB center in the absence of inter- 
band mixing may broaden into a narrow band of extended 
states when the effect of mixing of states of opposite chi- 
rality is taken into account. We also provide a physi- 
cal picture to show how the mixing of states of opposite 
chirality may possibly lead to the existence of extended 
state bands. Based on above results, we propose a new 
phase diagram in which a narrow metallic phase exists 
between two neighboring IQHE phases and between an 
IQHE phase and an insulating phase. This new phase di- 
agram is consistent with non-scaling behaviors observed 
in recent experiments. A standard scaling analysis on 
nonscaling experiment data2i supports our results. How- 
ever, due to finite-size effects, our numerical results can 
also be explained based on the assumption of a single 
critical point. Thus further study on large system size is 
still needed to conclude whether there are extended state 
bands in quantum Hall systems in thermodynamic limit. 



FIG. 19: A two-channel network model of 2 x 2 nodes with 
periodic boundaries along both directions. Z-s are the wave- 
function amplitudes on links. The notations are as follows. 
H and V stand for horizontal and vertical links, respectively. 
u (I) is for the upper (lower) LB. Si are SO(4) matrices de- 
scribing tunnelling at nodes, and Mi are U(2) matrices for 
interband mixing. 



In this appendix, we explicitly construct the evolution 
matrix U for a 2 x 2 two-channel CC-network model as 
shown in Fig ll9l Periodical boundary conditions in both 
directions are imposed as explained in section III. Z-s 
are the wavefunction amplitudes on links. The notations 
are as follows. H and V stand for horizontal and vertical 
links, respectively, u (I) is for the upper (lower) LB. 
Si are SO (4) matrices defined in Eq.Q describing the 
tunnelling on nodes, and Mi are U{2) matrices defined 
in Eqs. [S]and|S]describing interband mixing. From Fig ll9l 
we can obtain 
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where 1 and 6 are the 2x2 identity and zero matrices, 
respectively. If we define 
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The evolution operator U is 
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where is the 4x4 zero matrix. It has the structure of 
Eq.@. 
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